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We examine the spectral function of the single electron Green function at finite tempera- 
tures for the Tomonaga-Luttinger model which consists of the mutual interaction with only 
the forward scattering. The spectral weight, which is calculated as a function of the fre- 
quency with the fixed wave number, shows that several peaks originating in the excitation 
spectra of charge and spin fluctuations vary into a single peak by the increase of temperature. 

§1. Introduction 

The Tomonaga-Luttinger model§'''&, which can be solved exactly, is a funda- 
mental model for studying the electronic state in one-dimensional interacting electron 
systems. Since the excitations are gapless for both charge and spin degree of freedom, 
the response functions for the superconducting state, spin density wave state and 
charge density wave state show the power law behavior with resnect to frequency 
and temperature where the exponents depend on the interactions. eP'cP 

Noticeable properties of the Tomonaga-Luttinger model have been obtained in 
the state of the single electron. The power law dependence has been found in the 
momentujiLdistribution functionU^ around the Fermi momentum and the density of 
states B'' IIP 'S^ around the Fermi energy. These facts are characteristic of the Luttinger 
liquid & where the quantum fluctuation in one-dimension leads to the marginal be- 
havior between the metallic state and the ordered state. The detail of the electronic 
state is obtained in the spectral function which is calculated from the imaginary 
part of the single electron Green function. The spectral function at absolute zero 
temperature exhibits several peaks which are associated with the fluctuation of the 
pairing electron and the separation of the charge and spin degrees of freedom in 
the presence of the interaction. HSP'Q) The spectral function for the Hubbard model 
with the quarter-filled band and the infinite repulsive interactional) has been cal- 
culated where the frequency-dependence is similar to that of the Luttinger liquid. 
However the temperature-dependence of the spectral function of the Tomonaga- 
Luttinger model is not studied qualitatively although the thermal fiuctuation takes 
an important role for the electronic state around the Fermi energy. 

In the present paper, we examine the spectral function at finite temperatures to 
understand how several peaks are varied by the thermal fiuctuation. In §2, formu- 
lation is given in terms of the single electron Green function at finite temperatures. 
In §3, the frequency-dependence of the spectral function is examined numerically 
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with some choices of interaction and momentum. The crossover from the ground 
state to the state expected at finite temperature is demonstrated. §4 is devoted to 
discussion. 



§2. Formulation 

The Tomonaga-Luttinger model is given by 

r=±s=tl k 



r,s,s' ki,k2,q 



{92\\Ss,s' + 92±Ss-s') C^^^s^rCk2,s'-rCk2+q,s'-rCk:,- 



q,s,r 



(2-1) 



where C^^ ^ denotes the creation operator of the fermion with spin s =\ (i)(= +(— )) 
and the momentum k being positive (negative) for r = +(— ). The first term is the 
kinetic energy where vy and k^ are Fermi velocity and Fermi momentum respectively. 
In the second term, the quantity 52||,± (= fi'2||,±/'''"^F) denotes the normalized coupling 
constant of the interaction for the forward scattering between two kinds of electrons 
with r = + and r = — and ^4j_ (= (74_L/7rz;F) is that for the same kind of electrons. 
Quantities A~'^ and L denote the momentum cutoff of the interaction and the length 
of the system respectively. 

Based on the bosonization method, 11* eq. (2-1) in case of \q\ ^^1^^ is expressed 



as 



E ? 



v=p,(j 



1 



(2-2) 



Vu =vf\/(1± 54^/2) 



5211 ±92±)/2 



1 ± ff4±/2 - (g2|| ±g2±)/2 

^ 1 ± 54±/2 + (5211 ±52±)/2 ' 

where +(— ) corresponds to p{(t) and 



o,±{x) = ^E 



27r 



2^ Lp 
27r 



-e 2 



-iqx ^ 
k,s 



k+q,s,+^k,s,+ 



k+q,s, 



k,s. 



e 2 



\q\-tqx 



k+q,s,+^k,s,+ 



Ck,s,+ ± Cf,j^g ,, _Ck 



k,s 



(2-3) 
(2-4) 

(2-5) 
(2-6) 



Equations (2-5) and (2-6), which denote the phase variables for the charge and spin 
fluctuations respectively, satisfy the commutation relation given by 



(2-7) 
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From eqs. ( |2-2| ) and (2-7), one obtains the excitation spectra, Vpq and Vf^q for the 
charge fluctuation and the spin fluctuation respectively. By introducing the cut- 
off parameter ao(— ^ +0) for the convergence and making use of the fermion field 
operator defined asi^ 



exp 



(2-8) 



- [(0p,+(x, t) + r9p^_{x, t)) + s ie^,+{x, t) + r9^,_{x, t))] 

the retarded Green function at finite temperatures is calculated as& 
G^,(x,t) = -ie{t) (^'F^^,{x,t)^[^,{o,o)+ + 'F^^,{o,o)+^[^,{x,t)) ^ 



27r 

xS{x,t,T) + 



n 



A 



{ao + i{vyt - rx))^ \{A + ivytf + x"^ 



fx - 


^ —x\ 




.-t)_ 



■{x,t,T) 



oo 



n n 

v=p,cT n=l 



1 + 



' v^t — rx ~ 



1 



1 + 



^Vut — rx' 



1 + 



' Vi,t + rx\^ 



(2-9) 



-1 



(2-10) 



The quantity ji,, which denotes the magnitude of the interaction, is defined by 

7. = (r?. + r/;' - 2)/8 . (2-11) 

The quantity, T, is the temperature and is taken as unity. In deriving eq. (2-9) 
from eq. ( |2-8D , we have used an approximation that the second line of eq. (2-9) 
is correct in case of |x|,|z;Ft| ^» A and results in the extra factor of v-p / {ypVaY^"^ 
in case of |x|,|z;Ft| <C A compared with the exact one. Actually, by noting 
that S{x, t,T) — > 1 in the limit of absolute zero temperature, eq. (2-9) in case 
of T = becomes equal to the retarded Green function obtained by Luther and 
Pescheli^'S^'i' which is valid for the length scale being larger than A. Therefore the 
following calculation of the spectral function is justified when the frequency ( the 
momentum ) is smaller than vp/A {A~^). We note the difference between two kinds 
of cutoff parameters in the Tomonaga-Luttinger model. The quantity A leading to 
the momentum cutoff of the interaction plays an essential role for the existence of 
the characteristic energy, vp/A, but the quantity ao vanishes in the end by taking 
the limit of ao ~^ 0. 

In terms of eq. (2-9), the spectral function is calculated as 



Ar{q,uj) = Im 

TT 



POO POO 



dx I dte-*^(''+"'=F)^-^*>G,^(x,t) 

i{qx—LUt) 
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n 



1 



+ 



u=p,a V (ao + i{vyt - rx)) 2 V (yl + ivyt) + 



S{x,t,T) 



(2-12) 



Wc examine eq. (2-12) for two kinds of cases. The case (i) is given by (/2|| = 5 7^ 
and g2± = 94± = which corresponds to the spinless Tomonaga-Luttinger model. 
The case (ii) is given by g2\\ = g2± = 94± = 5 7^ which represents the spinful 
Tomonaga-Luttinger model. 

In the case (i), eq. (2-12) is rewritten as 



Ar{q,uj) 



Fi{f2i,T)F2{f22,T) + 



no 



-fix 



(2-13) 



By use of i7i = (cj + rvq) /2v, Q2 = iy — rvq) /2v, si = vt — rx and S2 = vt + rx, 
quantities Fi(fli,T) and F2{f22,T) are written as 



Fi{ni,T) 

-r- ' 



dsi . I - 

-00 CKo + IS\ \A + ISl 



A 



'y 00 

X 
n=l 



n 



1+ 



Sl 



Jo 



T 

Iv 



1 + 



Sl 



A+^ 



Sl (.12 + sl) 2 sinh {Tsi/v) 



1 + 



Sl 



A + ^ 



21 -T 



(2-14) 



F2{n2,T) 



J — ( 



ds2e 



A 



A + is2 



n 

n=l 



1 + 



S2 



A + ^ 



= 2A"< / ds2 



cos [l22S2-7tan-i(f)l g 



n=l 



1 + 



S2 



A + n^ 



{A^ + sir- 

where t> = 1;^ = = vy[1 - [92/2^]^'^ and 7 = [(1 - (52/2)^)-^/^ - l]/2. 
In the case (ii), eq. (2-12) is rewritten as 



(2-15) 



(27r)2|i;p - Va\ 

X n 

v=p,a 



dSa 



— CO 

A 



A 



+ 



{ao + isy)^ \A + isy) \A + i{aySp + Ksa) 



:{x,t,T) 



yl2(7p+7<,) 
{2Tr)'^\Vp - Va\ 



[Fs{np,f2^,T)+Fc{np,n^,T)] + 



-Qp 

-On 



(2-16) 
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By use of = {uj-rVaq)/{vp-Va), Qp = {rvpq-uj)/{vp-v„), Sp = Vpt-rx, = 
Vat - rx, ap = {vp + Va)/{vp - v^), bp = -2vp/{vp - v^), = 2va/{vp - v^), = 
— {vp + Va)/{vp — Va) , quantities Fs{fip, f2a,T) and Fc{f2p, f2a,T) are expressed as 



Fsinp,na,T) = 2 dsp dsa 



OD POO 



Yl s^{A^ + sly-- [A^ + {a^sp + b^Sa 



X sm 



v=p,a 



tan 



1 1 ^ 
A 



+ tan 



1 / CLpS p + bpSc 



A 



X ^{sp,Sa,T) , 



(2-17) 



Fc{f2p,f2a,T) = 2 ds 



dSa- 



Yl \s^\ (yl^ + sly-- {A'^ + {a^Sp + b^Sa] 



i/=p,a 



X COS 



^crSp -\- fIpS(j ^ ^ 



y=p,a 



tan 



I j 

~A 



+ tan 



1 f a-iySp + bySc 
A 



xC{sp,Sa,T) , (2-18) 

where 7p = [(1 - (25/(2 + 5))')"'/' - l]/4, 7<x = 0, z;^ = ?;f[(1 + g/2f - g^]^^ and 
t^o- = vf{1 — g/2l. In deriving eqs. (2-17) and (2-18), we made use of the following 



approximation, t3' 

six,t,T)^ n 



Tsp/vu 



V=p,(J 



sinh {Tsy/vu) 



i+7. 



TjauSp + bySa.)/vy 
sinh {T{aySp + bySa)/vy) 



^ 7i^ 



(2-19) 



which was obtained by discarding A in eq. (2-10). Such a treatment is valid for 
T^v-p/A in the present calculation. 

In addition to the case (ii) , we examine the case of 54^ 7^ and zero otherwise 
which corresponds to the one-branch Luttinger liquid, i.e., the forward scattering 
within the same kind of electrons. In this case, the spectral function A+(g,u;) is 
obtained by putting 7p = 70- = in eq. (2-16) where Aj^{q,uj) also shows the 
separation of the charge and spin degrees of freedom. 

§3. Spectral Function 

We evaluate the spectral weight of eq.(2-12) which is normalized as, 

A+{q,i:i) = A+{q,oj)vpA-^ . (3-1) 

Quantities g, u), v and T are also normalized as g = qA, uj = uj/{vfA~^), v = v/vp 
and T = T/{vpA~^) respectively. Prom eq.( p-ll ), the parameter for the interaction 
is defined as 



a = 2 (7p + 7a 



(3-2) 
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We examine A^{q,u)) in detail by choosing a= 0.125 which corresponds to the 
hmit of the large repulsive interaction for one-dimensional Hubbard model. 0^ Since 
Ar{q,uj) = Ar{—q, —23) = A^riQ, —23), we investigate numerically A^{q,u}) in case of 
q>0. 

3.1. spinless case 

The model in the case of g2\\ ^ and g2_L = gi± = is equivalent to that of the 
spinless fermion since there is no distinction between the charge fluctuation and the 
spin fluctuation, i.e., 7p = 7o-. 

First, A^(q,u!) with q = is examined. 
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1 1 1 1 1 

&A=0.125 


1 1 1 1 


tf=0 


-1^=0 
^0.05 


/ 
1. 


1^0.1 


'.^^^J^ 


VS<-.^0.2 

^-^i-~-V--ri 



-0.2 0.2 



Fig. 1. The tj-dependence of spectral function 
A^{ci,lS) of the spinless model where g = 0, 
a = 0.125 and f is chosen as f =0, 0.05, 0.1 
and 0.2. 



2 




1 2 



1 a 
T 

Fig. 2. The T-dependence of the imagi- 
nary part of the normalized self-energy, 
-lmr^(0,0) -ImX'^(0,0)yl/'UF) is 

shown in the case of a = 0.125 and 1.5. 
The dashed curve is the asymptotic value 
given by eq. (3-3) and the inset shows the 
coefficient given by eq. (3-4). 



In Fig. 1, u3) with a = 0.125 and g = is shown as a function of u) with 

some choices of T. In case of T = 0, the spectral function^^+(g, 23) divergesat 23 = 
where A^{q,uj) with q = and small |23| is given by 23) oc 23"~-'^.I13)'EP For 

T / 0, A+{q,uj) with 23 = becomes finite and the peak of A^{q,Lj) decreases and 
the width increases due to the thermal fluctuation. The half width of A^{q, 23), which 
is defined by Z\23, is proportional to T within the numerical accuracy of the present 
calculation, e.g. Auj ~ OAT (~ 2f ) for f < 0.1 in case of a = 0.125(= 0.5). The 
height of the peak at low temperatures is given by ^-i-(0, 0) oc T°^^ which is related 
to the imaginary part of the self-energy of the Fourier transform of the retarded 
Green function, eq. (2-9). Actually, by defining lmU^{q,uj) as the imaginary part 
of the Green function, one finds the relation that Imi7^(0, 0) = — [7r^(0, 0)]~^ where 
Rei:^(0,0) = 0. From eq. ( ^ ), the quantity ImZ'^(0,0) in case of T < vA'^ is 
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calculated as 



lmi;^(0,0) ~ C(a)f ^ 



-K^^^ sin (iT^) 



2v JO 

oo 



dy 



X / dy (sinh(7ry)) 
Jo 



sm 



h(vry)/ 



(3-3) 



(3-4) 



In Fig. 2, the normalized quantity of Imi7^(0,0) as a function of T is shown by 
the solid curves for a = 0.125 and 0.5. The dashed curve denotes the asymptotic 
value given by eq. (3-3) where the good coincidence between the solid curve and 
the dashed curve is obtained at low temperatures. The actual value of ImZ'^(0, 0) 
is smaller (larger ) than the asymptotic one value where the crossover takes place 
around a ~ 0.32 at low temperatures. In the inset, C(a) of eq. (3-4) is shown as 
a function of a. The quantity C{a) takes a maximum around a ~ 0.44 and C{a) 
becomes zero at a = 0, indicating the fact that Ar{q,Lo) — > 6{u! — vpq) at a — > 0. 

Next A^{^,uj) in case of q = 0.1 
is examined. The w-dependence of 
A^(q,ui) in case of q = 0.125 is shown 
in Fig. 3. The quantity A+{q,uj) at 
T = shows not only the main peak 
around u! = vq but also that near ui = 
—vq where A^{q,u)) oc {u! — 
for ui^vq, A^{q,uj) oc \u) + vq\°^^'^ for 



< 



- vq and A^{q,uj) = for \u!\ < vq. 
The peak around ui = —vq comes 
from the particle-hole excitations be- 

reen two kinds of electrons with r = it. 
The finite magnitude of Aj^{q,u;) ap- 
pears in the interval region of |a;| < vq 
at finite temperatures. By the increase 
of T, these two peaks are suppressed and 
merges into a single peak. 




Fig. 3. The (^-dependence of spectral function 
A+ {q, ui) of the spinless model where q = 
0.1, a = 0.125 and f is chosen as T=0, 
0.01, 0.02, 0.05 and 0.1. 



3.2. spinful case 

We examine A^{q,ui) in the presence of the interactions g2\\, g2± and g4± which 
result in the separation of the charge degree of freedom from the spin degree of 
freedom. When q = 0, the w-dependence of A^{q,uj) is similar to Fig. 1 since the 
excitation spectra of both charge fluctuation and spin fluctuation become equal to 
zero at q = 0. In Fig. 4, A^{q,Lo) in case of g = 0.5 and q;=0.125 is shown with 
some choices of T where A^{q,oj) with w < is multiplied by 10. In case of T = 0, 
there are several kinds of edges in A^{q,Ld) which originate in excitation spectra of 
the spin and charge fluctuations. Their asymptotic forms in the case of 7a- = are 
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fix 

fir 



0.2 



1 

&A=0.125 






- q=o.5 ,^ 

T=0.01 




> 


0.03- 






0.06- 
0.1- 










^ - 


J 






1 




1 



! A 
&X 



Fig. 4. The cj-dependence of spectral func- 
tion A+{q,ij) of the spinful model where 
q = 0.5, a = 0.125 and T is chosen as 
r= 0.01(solid curve), 0.03(dashed curve), 
0.06(dash-dotted curve) and 0.1 (dotted 
curve) . 



, , , 1 
&A=0.125 

T=0.03 


1 A 

L|— U. 1 




.1 . 0.3 


1 


ih;^ .-^0.5 
1 i; 1—' ■. 





1 

! A 
&X 



Fig. 5. The tj-dependence of spectral function 
A+ (q, ui) of the spinlful model where a = 
0.125, T = 0.03 and q is chosen as g=0.1, 
0.2, 0.3 and 0.5. 



given 



by0)|) 



(3-5) 



A+^a ~ ^+,p at T=0.03 and that 
dominant j4+^p into the dominant A. 



and A+{q,uj) = for — fpQ' < to < v^q ■ 

The T-dependence of A+{q,u}) is examined by defining A+^p (^+,0-) as tD) 
corresponding to the peak located near uj = Vpq {uj = Vo-q) where v^^q < Vpq. The 
result that A^^„ < ^+,p in case of T = 0.01 can be understood from the fact that, at 
T = 0, the exponent for the divergence of the charge excitation is larger than that 
of the spin excitation as is seen from eqs. (3-5). By the increase of T, one finds that 

> ^+,p for T=0.06. Such a crossover from the 
is characteristic of the finite temperature. 
The spin excitation, which has the energy lower than that of the charge excitation, 
has the large effect on A^{q,u!) and give rise to the dominant A^^p as is seen from 
the interference term, f2pSa + i^aSp, in eq. (2-16). The quantity A+{q,uj) at T = 0.1 
shows that the two peaks at low temperatures becomes a single peak with the broad 
width. When T increases, A^{q,ijj) with lij < decreases indicating the fact that 
the correlation by the interaction decreases by the thermal fluctuation. In Fig. 5, 
we show A^{q,ui) with the fixed T by choosing several q. The interval length of uj 
between two peaks, which exists for the large q, decreases by the decrease of q and 
vanishes for the small q, e.g., A^{q,Lij) for uj = 0.2 and 0.1. 
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For the comparison, A^{q,u!) with 
the interaction of only g4^± is shown in 
Fig. 6. In case of T = 0, A^{q,uj) 
exists only in the interval region of 
VaQ < u) < VpQ and show two singu- 
larities at the boundaries correspond- 
ing to the spin and charge excitations 
respectively. The exponents for the 
divergence at the charge excitation is 
the same as the spin excitation where 
Aj^{q,ui) (X {ui — Vaq)~^^'^ for 



{—uj + Vpq) for 



and A^{qS)) = for a; < v^q and 
UJ > Vpq .u) The effect of the thermal 
fluctuation, which leads to the crossover 
from two peaks to the single peak, is 
similar to Fig. 5. When uj = v^q or 
UJ = Vf^q, eq. (2-12) is rewritten as 



A. 



\q,v„q) = 



A+{q,Vpq) 



T 



-•La 



^=1 




^0 

1 — u 


: j 


\ I 


0.05 

—0.1 
\— 0.2 

v.. 


///' 

r-:U-^', 









&^ 


1 



Fig. 6. The (^-dependence of spectral function 
A+{q,uj) with the interaction, g4±, where 
q = 0.5, (?4 = g4±/{T^VF) and T is chosen 
as f =0, 0.05, 0.1 and 0.2. 



dy- 



cos 



dy 



VpVaTT 



where Ci = ^ 
and eqs. (3-7) oc T' 
rapidly compared with A_^ 



(sinhy) 2 

cos {^y) 
1 

(sinhy) 2 

1 

' 2 



(3-6) 
(3-7) 

'1/2,,- 



-1/2, 



Va\ ""^ /g°° dy (sinh y) 2 . Since eqs. (3-6) oc T' 
at low temperatures, it turns out that A^^p decreases 



§4. Discussion 



We examined the spectral function, A^(q,uj), of the Tomonaga-Luttinger model 
at finite temperatures for both the spinless case and the spinful case by choosing 
a = 0.125 which corresponds to the large limit of the repulsive interaction of the 
Hubbard model. In case of g = 0, A+{q,uj) shows the peak around uj = 0. By the 
increase of temperature, the height decreases and the width increases due to the 
thermal fluctuation. At low temperatures, A^(q,uj) with q = and uj = 0, shows 
the power law expressed as A^{q,Ld) oc T°'~^. In case of g > 0, there are two peaks 
located in the region of u; > and that of u; < for both the spinless and spinful 
cases. With increasing temperature, these two peaks moves to a single peak around 
uj = vq. In the spinful case, the peak with g > is separated into two peaks where 
the peak with the large uj, ^+,p, corresponds to the charge fluctuation and the peak 
with the small uj, ^+,0-, corresponds to the spin fluctuation. In the limit of low 
temperatures, one obtains ^+,0- < ^+,p, while one finds ^+,0- > ^+,0- by the increase 
of temperature. The latter result, which is explained by eq. (2T6), can be also 
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understood in terms of the general formula0) 
Mk,^) = y E [l(C^fc'')n,m|'e-*(^™-^^'")(l + e-T)S{En - - fi - co)] , (4-1) 



m,n 



where Z = J2n e^^"~^^"'^^'^ and Nn is the electron number. Quantities En and /x are 
the energy of the n-th eigenstate and the chemical potential respectively and {C'^)n,m 
denotes the matrix element between the n-th eigenstate and m-th eigenstate. Since 
VpQ > VcrQ in the present case, the factor, (1 + e~^/'^), with oj = v^q is larger than 
the factor with u; = Vpq. Therefore the peak for the spin fluctuation becomes larger 
than that for the charge fluctuation. At higher temperatures, these two peaks also 
merge into a single peak. 

Finally, we comment on the experiment on K0.3M0O3 where the angle-resolved 
photoemission spectroscopy reveals the two peaks corresponding to the charge and 
spin separations. HIP As for two peaks indicating the properties of the Tomonaga- 
Luttinger liquid, the peak with the lower energy is larger than the peak with the 
higher energy. In the Tomonaga-Luttinger model, the peak with lower energy is 
rather suppressed when the interaction for the charge density is strong enough. Then 
VoitllS' claimed that the peak with the higher energy is rather suppressed by the 
backward scattering. 113) Here we comment another possibility that the suppression 
of the peak with higher energy is attributable to the effect of thermal fluctuation as 
is found in Figs. 4 and 6. 
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